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Abstract. In this paper, we give an upper bound for a probabilistic distance 
between a Gaussian vector and a vector of U-statistics of Poisson point pro- 
cesses by applying Malliavin-Stein inequality on the Poisson space. 



Contents 

1. Introduction 1 

2. Malliavin Calculus on the Poisson space 2 

3. Multi-dimensional Malliavin-Stein inequality 4 

4. Application for multi-dimensional U-statistics 7 
References [n 



1. Introduction 

The theory of Malliavin calculus on the Poisson space was firstly studied by 
Nualart and Vives in their excellent paper of Strasbourg's seminars (NV90| . For 
some important contributions and applications for Poisson point processes, we refer 
to [HPA95t IWuOO[ IPTllI iLPTl] . The combination of Stein's method and Malli- 
avin calculus on the Poisson space related to the normal approximations of Poisson 
functionals has been considered by Peccati, Sole, Taqqu, Utzet and Zheng in their 
recent papers [PSTUlOi IPZTO] . 

The basic theory of U-statistics was introduced by Hoeffding |Hoe48j as a class 
of statistics that is especially important in estimation theory. Further applications 
are widely regarded in theory of random graphs, spatial statistics, theory of com- 
munication and stochastic geometry, see e.g. |Lee90[ IKB941 lBor96| . 



Recently, the idea of central limit theorems for U-statistics of Poisson point 
processes by using Malliavin calculus and Stein's method has been given by Reitzner 
and Schulte in |RS11| . Our main work in the current paper is to extend their 
results to vectors of U-statistics by applying the muti-dimensional Malliavin-Stein 
inequality, that was proved by Peccati and Zheng in |PZ10] . Some preliminaries 
of Malliavin calculus on the Poisson space and U-statistic will be introduced in 
Section 2. An upper bound for a probabilistic distance between a Gaussian vector 
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and a square integrable random variable with finite Wiener-Ito chaos expansions 
wiU be shown in Section 3 and its apphcation for multi-dimensional U-statistics of 
Poisson point processes will be given in Section 4. 

2. Malliavin Calculus on the Poisson space 

Let {E, A, jj) be some measure space with cr-finite measure /x. The Poisson point 

process or Poisson measure with intensity measure /i is a family of random variables 
{N{A)}AeA defined on some probability space (fi, J^, P) such that 

1. VA e A, N{A) is a Poisson random variable with rate ijl{A). 

2. If sets Ai, A2, . . . , An e A don't intersect then the corresponding random 
variables N{.) from i) are mutually independent. 

3. Vw e fi, N{.,uj) is a measure on {E,A). 

Note that, for some a-finite measure space {E,A,n), we can always set 

n 

n = {oj = ^6;,^, n G N U {00}, Zj e E}, 

where Sz denotes the Dirac mass at z, and for each A G A, we give the mapping A'' 
such that 

uj ^ N{A,iu) = iu{A). 

Moreover, the a-field J- is supposed to be the P-completion of the a-field generated 
by AT. 

Let A4{E) denote the space of all integer- valued ct- finite measures on E, which 
can be equipped with the smallest a-algebra E such that for each A E A then 
the mapping ij E Ai{E) 1— > ri{A) is measurable. We can give on the 
probability measure Pn induced by the Poisson measure N and denote LP{Pn) as 
the set of all measurable functions F : M{E) — )• M such that E[|F|p] < 00, where 
the expectation takes w.r.t. the probability measure Pn- 

Let LP{fjb^) be the space of all measurable functions / : — )• M such that 

11/11= j |/(zi,...,^n)|V"(c«^l,---,t^^n) <0O. 

Note that Z/^(/i") becomes a Hilbert space when we define on it the scalar product 

{f,g) = j fizi,...,Zn)g{zi,...,Zn)l^'\dzi,...,dZn). 

We denote Lfyjjj(/x") as the subset of symmetric functions in in the sense 

that the functions arc invariant under all permutations of the arguments. Now, 
for each A E A, we define the random variable N{A) = N{A) — ij,{A), which 
is also known as a compensated Poisson measure. For each symmetric function 
/ e Lgyjjj(/x"'), one can define the multiple Wiener-Ito integral /«(/) w.r.t the 

compensated Poisson measure A'' denoted by 

In{f)= [ f{zi,...,Zn)N''{dZu...,dZn). (2.1) 
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At first, let 5„ be the class of simple functions /, which takes the form 

m 

f{zi,Z2,--,Zn) = ^>'kl^m ^^^^^^(k){zi, Z2, Zn), (2.2) 
k=l 

where Af^^ G .4, Afc € M and the sets A[''^ x . . . x yll*^'' are pairwise disjoint such that 
/ is symmetric and vanishes on diagonals, that means /(^i, . . . , Zn) = if = Zj 
for some i ^ j. The multiple Wiener- Ito integral for a simple function / in the 
form (|2.2p with respect to the compensated Poisson measure N is defined by 

m 

^«(/) -I]AfciV(Ai)...iV(A„). 

Since the class 5„ is dense in L'^y^{fi") then for every / e Lgyj^(^"), there exists 
a sequence {fi}i>o C 5„ such that fi ^ f in Lgy,,^(^"). Moreover, one can show 
that E[/„(/;)^] = Hence, the the multiple Wiener-Ito integral /„(/) for 

a symmetric function / e L^y^(^") can be defined as the limit of the sequence 
{In{fi)}i>o in L'^{Pn) and we denote it as (j2.ip . 

Proposition 2.1. For n,m> 1 and f G L'^y^{fi"), g G Lgyj^(^™), i/ien 

(1) E[/„(/)] -0, 

(2) E[/„(/)/„(.g)] =<5 

where dn,m is the Kronecker delta. 

For a measurable function F : A^(£') — > M and z E E we define the difference 
operator as 

D,F{r^)^F{r^ + S,)-F{rj), 

where dz is the Dirac measure at the point z. The iterated difference operator is 
given by 

Dz,,...,z„F = Dz,Dz,^,„^z^F. 
We define the kernels of F as functions /„ : E" — ^ R given by 

/„(zi,...,z„) = ^E[Dzi„,,^Zr,F],n > 1, 

Note that /„ is a symmetric function. 

We define the Ornstein-Uhlenbeck generator as 

LF{r^) = J (F(77 - Sz) - F{Tj)Mdz) - J {F{f^) - F{tj + S,)) fi{dz). 

E E 

Proposition 2.2. For each F e L'^{Pn), then the kernels /„ are elements of 
n > 1 and uniquely admit the Wiener-Ito chaos expansion in the form 

oo 

F = E[F]+J2ln{fn), 
n=l 

where the sum converges in L'^{Pn). Furthermore, for F,G £ L'^{Pn) 

oo 

Cov(^^,G) = E[FG] - E[^^]E[G] = «!(/„, 5„)l^(m")- 
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Proposition 2.3. Let F e L^{Pn), and assume that 

oo 

J2nn\\\fnf < oo. 

n=l 

Then the difference of F at z € E is given by 

oo 

D,F=J2nIn-l{fn{y,-))- 

n=l 

Proposition 2.4. For each random variable F E L?'{Pm) such that 

oo 

^z2z!||/„||2<oo, 

n=l 

then the Ornstein- Uhlenbeck generator L is calculated as 

oo 
n=l 

Moreover, its inverse operator is calculated as 

oo 

L-'F=-Y,-lMn). 
n—l 

for each F e L'^{Pn) such that E[i^] = 0. 

For more details of the Malliavin calculus on Poisson space, we refer the reader 
to [NV90] . 

3. Multi-dimensional Malliavin-Stein inequality 

In the next sequence, we use the probabilistic distance of two d-dimensional 
random vectors X,Y such that E(||X||jjd), E(||y||ud) < oo, which is defined by 

A(X,y) = sup|E(g(X))-E(5(r))|, 

gen 

where Ti, is the family of all real- valued functions g €E C^(R'^) such that 

bllLip = sup '^["^^'-f^^' < 1, sup ||Hess(ff(x))|l < 1. 
In the above inequality, 

stands for the Hessian matrix of g evaluated at a point z and we use the notation 
of operator norm for a d x d real matrix A given by 

\\M = sup IjAxllRd. 

Nll3id=l 
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Theorem 3.1 (Multi-dimensional Malliavin-Stein inequality). Consider a random 
vector F (Fi, . . . , Fd) C L'^{PN),d > 2 such that for I < i < d, Ft e dom{D), 
and E(_Fi) = 0. Suppose that X ^ Afd{0, C), where C = {C(z, j) : i, = 1, . . . , c?} is 
a d X d positive definite symmetric matrix. Then, 



A{F,X)<\\C-^\\Ct^\ 



^ n{C{i,j) - {DF-DL-^F,)L.^^)y 



-\\C-'\nC\\ / Kdz)E 



Y^\d,fA [Y.\d,l-'fA 

\i=l / \i=l ) 



For the proof, we refer to |PZ10) . 



Now we consider a d-dimensional random vector F = (F\, . . . ,Fd) C L'^{Pn) 
with the covariance matrix E = j) : i, j — 1, . . . ,d} such that each component 

(n) 

Fi has finite Wiener- Ito chaos expansions with kernels , which vanishes if n > fc. 
Let give the centered random vector 



G = \/CE-i(F-E[F]), 

where VA stands for the square root of a positive definite matrix A, i.e if A has 
the eigenvalues decomposition A = P^^diag(Ai, . . . , Xd)P then 

VA = P-Miag(v/A7, . . . , ^/x'd)P. 
Let us use vector notations 

V,F^{D,Fu...,D,Fd), V,LF={D,LFi,...,D,LFd) 
and note that the inequality 

trace(AB) < trace(A)trace(B) (3.1) 
holds for all positive definite matrices A, B. 

Therefore, by the properties matrix trace trace(yli?) ~ trace(i?A) and using the 
inequality (|3.ip . we have 

2' 



— trace 



trace 



C- J^fiidz)VM~V,LGf^ 
VCS-i " / m(c^^)V^[F - E(F)](-V^L[F - E{F)]f^ 
{C^-'f (e - /i(dz)V,[F - E(P)](-V.L[P - HfW) 
E - / M(dz)V,[F~E(F)](~V,L[P-E(P)]) 

J E 



< trace[(CS-i)2]trace 



where 



A\\f = A/tracc(A^A) 
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denotes the Frobenius norm of matrix A. 
Note that 



S(i,j) = Cov(i^,,F,) = ^n!(/f\/j") 



n=l 



and 



{D[F,-^{F,)],-DL-'[Fj-^{Fj)])r..^,^ = ( £ nln-,{ft\z, •)), In-i{ff{z, ■)) 



\n=l 



Hence, 



/ k 



< 



k 



n=l 



n!(/f\/i"^)^,(^„^-n(/n-l(/i"^(^,-)),/n-l(/|"^(.,-)) 



1,2 (^) 



+ ^ E n-(l„.,{fl-\z,-))J^.,{f^"'\z,-))) 

n,m=l,n7=m 

= fc^ 5: n2Var('(/„_i(/(")(.,.)),/™-i(/f^(.,-)))^,^J. 



l<n,m<fe 

It follows that 



^ E[(C(z,j) - (i?G,,-Z?L-iG,)^2(^))2] 
ll^^-'II^A E E Var((j„_i(/(")(.,.)),/™-i(/f)(.,-)))^ ). 



(3.2) 



Morever, by using Holder inequality and the property of matrix norm, we have 



/Lt(dz)E 



) \i=l / 



<d3/2 / Mdz)E[||V.G||2,||V.L-iG||i 



= d^l'' [ M(rf^)Ef||VGS^V.[F-E(F)]||2,||VGS^V.L-i[^-E(F)]|U. 
<rf=^/2||^/GS^f / Mrf-^)E[||V.[F-E(F)]||^.||V.L-Mi=^-E(F)]|U.] 

j2mF,-E{F^]A 



1 ^ 



X / ij,{dz)E 



J2\D.L-'[F,-E{F,)f 



1/2 



< 



1/2 



MALLIAVIN-STEIN METHOD FOR MULTI-DIMENSIONAL POISSON U-STATISTICS 7 



^ (J2 /^M(d^)E [\D,L-'[F,-I)iF,)]\']j 
<d^Vc^r (E / /i(dz)fc3^n4E[/„_i(/(")(z,.))^ 

n=l 



1/2 



1/2 



n=l J 



1/2 



< rf2fc^/2||ycS^||3(trace(E)) 



d k 



'^\EEE[ii^-i(/f^(-.-))^ip' 

\ i=l ri=l 



(3.3) 



Substituting p. 21) and p.3p to the inequality in Theorem 3.1 for G, we obtain that 

Theorem 3.2. Let give a d- dimensional Gaussian random variable X ^ Md{0, C). 
Assume that F = (Fi, . . . ,Fd) C L'^{Pn) such that Cov(Fi,Fj) = 'S{i,j), i,j = l,d 
and Fi has finite Wiener-Ito chaos expansions with kernels /|"'' , which vanishes if 
n > k. Then 



A (VCS-i {F - F,{F)) , x) < 



2tt 



+k^C^-^F\\C~^\\C\ 



1/2 



3/2||C||(trace(E))i/2^ 


EEe[||/„_i(/(")(z,.))^IP^ 

i—l n—1 




d k 

E E Var 

z,j — 1 n,m— 1 


^(/„_i(/(")(z,.)),/«-i(/rH-,-)))^,^^j 



4. Application for multi-dimensional U-statistics 

In this section we consider the d-dimensional vector of U-statistics of the Poisson 
point process N 



F : 



E (;ii(zi,...,zfcj,..., E (l^dizi, . ■ . , ZkJ \ , 

^{zi,...,Zk^)eSki{N) {zi,...,ZkJeSk^{N) 



(4.1) 

where 0i G Lly^{^^'), and ^fc. (TV) denotes the set of ah fc^-tuples of distinct points 
of N. This means that each component 

Fi = E (Piizi, ■ . ■,Zk,) 

(zi,...,2^)GSfc,(A') 

is an U-statistic of order ki with respect to the Poisson point process N, i ^ l,d. 

The fohowing properties of (one-dimensional) U-statistics are obtained by Re- 
itzner and Schulte in [RSllj 



8 



NGUYEN TUAN MINH 



Proposition 4.1. Let F G L'^{Pn) be a U -statistic of order k in the form 

F= ^ 0(^1, . . . ,Zfe) 

(zi,...,Zfc)6Sfc(7V) 

Then the kernels of the Wiener-Ito chaos expansion of F have the form 
'k'^ 



/„(zi, . . . ,z„) 



(f>[zi,...,Zn,Xi,...,Xk-n)tJ''' " (dxi , . . . , dXfe-n ) , n<k 



n> k. 



Proposition 4.2. Assume F G L'^{Pisi), then 

(1) If F is a U-statistic, then F has a finite Wiener-Ito chaos expansion with 
kernels /„ € n n = 1, . . . 

(2) If F has a finite Wiener-Ito chaos expansion with kernels fn G H 

n ~ 1, . . . ,k, then F is a finite sum of U-statistics and a constant. 



Proposition 4.3. Let fi e S^i, i — 1, . . . ,m and 11 be the set of all partitions of 
Z = {z[^\ . . . , Zni , ■ • ■ , z[™\ . . . , zii^}, Hi < ki such that for each tt € 11, 

(1) Zi'^\ Zf^"' Cz Z , I h are always in different subsets o/tt, and such that 

(2) every subset of tt has at least two elements. 

For every partition tt S 11 we define an operator that replaces all elements of Z 
in YVl^i fi{'^i \ ■ • ■ 7 Zn}) that belong to the same subset of tt by a new variable Xj, 
J = 1, . . . , |7r|, where \'k\ denotes the number of subsets of the partition tt. Then 



E 



/■nri 
R''{\\_fi{-)){xi,-- . ,a;|^i)//l''l(da;i, . . .,dx\^\). 



Using the Proposition 4.3 and the same technique in |RS11) (Lemma 4.6), we 
also obtain that \i F ^ {Fi, F2, ...,Fd) C L^{Pn) is a vector of U-statistics in the 
form (|4.ip such that (f)i,i = l,d are simple functions, then all kernels /f""* are also 
simple functions and 

Var(^(/„_i(/(")(.,.)),/™-i(/j"'(z,-)))^,^^J 



E 



(l,.^,{ft\z,-)),Ir.Mf^^"'\z,-)))[^^^^ 



E 



/„_i(/(")(z,.))%.-i(/r(y,-))' t^''idy,dz) 



and 



(|/.'"^(0/f^(0/r(-)/r(-)|)(^i,.-.,^MW'^i(d^i,-.-,rf^i.i), 



E 



\\in-i{ft\z,-)r\ 



E 



R 



^(|/f^(.)/i"^(-)/i"^(-)/i"'(-)|)(xi,.-.,^MVl^l(d^i,..-,dxM), 
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where Tln,n stands for the set of partitions satisfying the conditions in Proposition 

zt-^wHhZ -i^^^ 7(1) 7(2) (2) (3) (3) (4) (4) l.„^jfT ^ 

^n,m denotes the set of all partitions in Iln,m of such that for any tt e Hn.m and any 
decomposition of {1, 2, 3, 4} into two disjoint sets Mi, M2 there are i E Mi,j € M2 
and two variables zj^ , z," which are in the same subset of tt. 
By the formula of kernels in Proposition 4.1., we note that 



(yi,---,2/|7r|), 



^K|+2(;c.+;c,-,-.0(rf^(i)^ . . . , d4'L' '^^f ^ ■ • ■ dxi%^,dy,, ...dy^^^). (4.2) 
This fact follows that 

Theorem 4.1. Assume that F = {Fi, F2, Fd) C L'^{Pn) is a vector of U- 
statistics in the form CT71) such that — are simple functions. Then 



A (VCS-i {F - E(F)) , < 
^d^fc^/^llVcS^f ||C-i||3/2||C||(trace(E))i/^ 



+e\\cj:-'y\\c-'\\\\c\\'/\ 



d k 



\ EE^^"'»(^'*) 

\ i=l n=l 



d k 



A E E ^n,m{i,j), 
\ n,m— 1 



where k = max{fci,l < i <,d} and M„^m(i,j'),l < i,j < d,l < n,m < k are 
defined in (4^- 



Now, we consider that F — (i^i, F2, Fd) C L^{Pn) a vector of {/-statistics in 
the form l|4.ip such that 

E \Mzi,...,zk,)\eL\PN). 

(zi,...,Zfc.)eSfc^(Af) 

Then, for each i = l,2,...,d there exists a sequence {(j)i.i}i>o C 5^. such that 
< and converges to (pi /I'^'-almost everywhere. Let give the vector of 
U-statistics i^*-'' — (i^i,;, . . . ,Fd,i), where 

FiJ = E 4'i.i{zi,- ■ ■ , Zki). 

(zi,...,ZfcjGSfc,(Ar) 

Hence, 

\F^,l\< E \'I^.A^l,---,Zk,)\< E |</'^(^l,...,^fcJ| ei'(Piv). 

(zi,...,2fc.)eSfc.(Af) (zi,...,ZfcjeSfc.(JV) 
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Its follow that Fi^i e L^{Pn), Fi^i converges to Fi almost surely and all kernels f^^"^ 
in the Wiener-Ito chaos expansion of Fi^i are simple functions. Note that 

nhi)^GoY{F,,Fj) = 

= (^^j) j J 4>t{zi-,----,Zn,Xi,...,Xk,-n)l^'^'~'^{dxi,...dXk^^n) 

(f>j{zi, ...,Zn,Xi,.. . . .dXfc^_„)^"((izi, . . .,dZn). 

Moreover, the integrals 



{(piizi, . . . , Z„,Xl, . . . ,Xk,-n)\fJ-''' "(da;i, . . . dXki-n) 



I m\^l, . ■ . ,^n,Xl, . . . ,.Lk,-n)\l-'^" 



\<j)j{zi, ...,Zn,Xi,.. ■,XkJ-n)\^J■ ' "{dxi, . . . da;fc^ _„) ^t" (dzi , . . .,dzn) 

always exist for 1 < n < ki,l < i,j < d. Therefore, by applying the Lebesgue 
dominated convergence theorem, we obtain that and E(i^i,;) 

E(i^i) for I -> oo. Hence, 

y^C(S('))-i (i^(') - E(F('))) ^ VCS-i [F - E(F)) 

almost surely for Z — ?> oo. Note that, the almost sure convergence implies the 
convergence in the probabilistic distance A and \Mn}m{i,j)\ < \Mn.m{hj)\ , where 
i) is defined when we replace 4>i,4>j by (j)i\4>^^ in (|4?2l) . Therefore, by 
using Theorem 4.1 and applying the triangular inequality, we conclude that 

Theorem 4.2. Assume that F = iFi,...,Fd) C L^{Pn) is a vector of U-statistics 
in the form J^.ip such that 

\Mzu---,zk,)\eL''{PN). 

izi,...,Zk,)eSk,{N) 

Then 



VA;^/2||%/CE^||3||C-i||3/2|lC|l(tracc(E))i/2^ 



d k 



\ i=l ,1=1 



+fc2|lCE-i|l^|lC-i|l|lC|li/2 ^ ^ M„,„,(z,j), 

\ ij—l n,m—l 

where k — max{fci,l < i <,d} and Af„^m(i,i),l < i,J l£ d,l < n,m < k are 
defined in (RTll). 



Corollary 4.3. Assume that {i^''''}/>o is a sequence of vectors of U-statistics, 
which are defined as in Theorem 4-2, such that 



max 
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for I -> oo, then the law of y/CiTWy^ (F^) - E(i^('))) converges to the multivari- 
ate Gaussian law AfdiO^C)- 
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